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We consider the production of charmed baryons and mesons in the proton- 
antiproton binary reactions at the energies of the future PANDA experiment. 
To describe these processes in terms of hadronic interaction models, one 
needs strong couphngs of the initial nucleons with the intermediate and final 
charmed hadrons. Similar couplings enter the models of binary reactions 
with strange hadrons. For both charmed and strange hadrons we employ 
the strong couplings and their ratios calculated from QCD light-cone sum 
rules. In this method finite masses of c and s quarks are taken into account. 
Employing the Kaidalov's quark-gluon string model with Regge poles and 
adjusting the normalization of the amplitudes in this model to the calcu- 
lated strong couplings, we estimate the production cross section of charmed 
hadrons. For pp — )• AcAc it can reach several tens of nb at piab = 15 GeV, 
whereas the cross sections of Sc and D pair production are predicted to be 
smaller. 



1 Introduction 



There is a vivid interest in the cross-section of charmed hadron production in the proton- 
antiproton colhsions to be measured by the future PANDA experiment (see, e.g., [1]). 
The amount of produced charmed mesons and baryons is important for assessing the 
abihty of this experiment to perform flavour-physics oriented studies, such as the mea- 
surement of charm-anticharm mixing, the search for CP-violation in D decays or the 
studies of Ac decays. A reliable estimate of the pp — )■ charm cross section is how- 
ever a very difficult task. The main problem is that the projected energy range (with 
the cm. energy ^/s varying from 2.25 GeV to 5.47 GeV), being reasonably high for 
the proton-antiproton collisions, is still not far from the threshold of charm-anticharm 
production. Several models of charm production at these energies can be found in 
the literature [2l[3ll3l[5l[6l[71[8l[9l [10] , their predictions differing by several orders of 
magnitude, as emphasized, e.g. in [8]. Especially difficult is to predict the inclusive 
charm-anticharm cross section in the situation where not many exclusive channels are 
open. Hence, it is more realistic to assess the exclusive production of baryons or mesons, 
such as pp — 7- AcAc, DD. A successful model of strange-hadron pair production in pp 
collisions, which was measured at similar energies, can serve as a useful tool, provided 
there is a reliable way to replace the model parameters of strange hadrons by the ones 
for charmed hadrons. The key parameters in many hadronic models of these processes 
are the strong couplings of strange or charmed baryons with mesons and nucleons. To 
relate them, the SU{4) j^-symmetry is frequently used in the literature. Note however, 
that it is difficult to justify this symmetry in QCD, due to the large mass difference of 
the c- and s- quarks, nic — mg » Aqcd- 

In this paper we employ the strong baryon-meson couplings of charmed and strange 
hadrons calculated from QCD light-cone sum rules (LCSR), where finite masses of c and 
s quark are taken into account. Recently, we calculated the charmed baryon strong 
couplings with a charmed meson and a nucleon. In addition to these results, here we ob- 
tain the corresponding strong couplings of strange hadrons. The nonperturbative inputs 
in the LCSR method are the universal nucleon distribution amplitudes (DA's). Hence, 
the extension of our calculation from charmed to strange hadrons is straightforward 
and is reduced to a replacement of the virtual c-quark by an s-quark in the underlying 
correlation functions. In what follows, we also employ the ratios of calculated strong 
couplings which are predicted from LCSR with smaller uncertainties than the individual 
couplings. 

The results for the strong couplings presented here can be used in various models of 
exclusive charm and strange hadron production. As an application of our results, we 
use the quark-gluon string (QGS) model of binary reactions developed by Kaidalov and 
his collaborators \12 \ \1?> \ [T ^ 115). One version of this model was already applied in [2J to 
estimate the charm production cross section in proton-antiproton collisions. We refine 
this model by introducing the helicity amplitudes and adjusting the two independent 
strong couplings to the LCSR estimates. 

In what follows, in Sect. 2 we present the LCSR results for strong couplings of charmed 
and strange hadrons. In Sect. 3 we demonstrate how the QGS model works for relatively 
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simple processes of meson pair production and trace the relation between the model 
parameters and strong couplings in QCD. In Sect. 4 we use the QGS model for pp 
binary reactions with charmed and strange hadrons, employing the strong couplings from 
LCSR and predict the charm production cross sections. Sect. 5 contains the concluding 
discussion. The two appendices contain: (A) the formulae for helicity amplitudes and 
(B) the derivation of the absorption factor in pp binary processes. 

2 Strong couplings from QCD light-cone sum rules 

The strong couplings of the Ac-baryon with the nucleon and D or D* meson are formally 
defined as the following hadronic matrix elements: 

{A,{P-q)\D{-q)N{P)) = gA^NDUA^P - q)il5UN{P), 
(A,(P - q)\D*{-q)N{P)) = ua^P " <?) (^91nd4 + i ^^^^ ^A^'-^^'^) ^n{P). (1) 

Note that the above couplings are defined in [IT] as residues at the D^*^ and Ac poles 
in double dispersion relations for the correlation functions with on-shell nucleon state, 
hence all three hadrons are on their mass-shell. The same definitions are valid for the 
Sc-baryon couplings as well as for the corresponding strange hadrons with the following 
replacements: Ac(Sc) — )• A(I1) and — )• iC^*) in the above. 

The AcND^*^ and S^iVD^*) strong couplings were calculated from LCSR in [11], where 
one can find the detailed description of the sum rule derivation. The results for the strong 
couplings which will be used in this paper are collected in Table [TJ Note that in p.lj two 
different interpolating currents for Ac and Ec baryons were used. With the procedure 
of eliminating the negative parity baryons suggested in that paper, the results agree 
within the uncertainties. In this paper we will only use the strong couplings obtained 
for the pseudoscalar interpolating current for A(c) and loffe current for S^c), respectively, 
because the sum rules in these cases have a comparatively lower background of higher 
states. In Tabled] also the ratios of strong couplings obtained from LCSR are presented, 
generally they have smaller estimated uncertainties, because of the common inputs used 
in the sum rules. 

Turning to strange hadrons, we employ the same LCSR method as in |161 [T7] and, 
replacing c-quark with the s-quark in the correlation function, calculate the AA^i^^*) 
and TiNK* couplings. The inputs used in LCSR consist of universal nucleon DA's 
which are taken from jl8j and explained in detail in [llj . In particular we use for the 
virtual c quark in the correlation function the value mc{mc) = 1.28 it 0.03 GeV. The 
flavour-specific input parameters which we adopt here for the sum rules involving strange 
hadrons are: the strange quark mass ms(2 GeV) = 98 it 16 MeV and the renormalization 
scale = 1-0 ± 0.2 GeV. Furthermore, one and the same range = 2.0 zt 0.5 GeV^ 
of the Borel parameter in the E and A channels is adopted, whereas for the K* channel 
we use M"^ = 1.0 it 0.5 GeV^. The threshold parameter in the LCSR for A(S) strong 
couplings is taken sq = 2.55 ± 0.10 GeV^ (sq = 2.75 ± 0.10 GeV^). The criteria of 
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Table 1: Numerical results for the strong couplings of charmed 11 If and strange baryons 
and their ratios obtained from LCSR with nucleon DA's. 

choosing the input parameters and the quark-hadron duahty ansatz in LCSR are the 
same as the ones used and discussed in ^llj. The two-point QCD sum rules for the 
decay constants of A and S baryons -with pseudoscalar and loffe currents respectively, 
are taken from [19]. Using the same definitions and notation as for the decay constants 
of charmed baryons in |11]). we obtain: 

Af ) = (0.87;°;?3) X 10-2 GeV^ , X^^ = x 10"^ GeV^ . (2) 

The resulting estimates of the strange baryon strong couplings and their ratios obtained 
from LCSR are presented in Tabled) Note that Ac and Sc belong to different SU (3) mul- 
tiplets (as opposed to A and S), and this circumstance explains a substantial difference 
between the ratios of tensor and vector strong couplings for these baryons. 

In the literature (see e.g., [3]), the strong couplings of charmed hadrons are estimated 
assuming SU (4) fi symmetry and equating dimensionless couplings, for example, assum- 

V(T) V(T) 

ing Vd* — 9kNK*- Such symmetry relations are difficult to justify from the point of 
view of QCD. Indeed, because of the large mass difference of c and s quarks, the kinemat- 
ical factors: masses and four-momenta entering the complete hadronic matrix elements 
differ significantly. In our approach we are not relying on any form of the SU{^)fi 
symmetry. Still, it is interesting to compare the strong couplings for the charmed and 
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Figure 1: The planar diagram of charmed baryon (a) and meson (b) pair production in 
pp collisions. 

strange baryons obtained from LCSR and collected in Table [TJ We find that the values 
of the dimensionless couplings are in the same ballpark, whereas there is a signif- 
icant difference between gj^jqj^* and ^JcAf-D*- '^^^ strange baryon couplings were also 
calculated with the Nijmegen potential model |20j of low-energy scattering, assuming 
SU (3) j^-symmetry. Expressed in terms of the dimensionless (7-couplings defined in ([T]) 
the results of [20J with their sign conventions are: 

= 13.4 - 17.5 , =-(4.3-6.1), ^liVi^* = 12.4 - 16.3, 

5EJVi^ = -(4.1-5.3), 5^^;^* = -(2.4-3.5), ^siVK* = -(1-3 ^ 4.6) . (3) 

Comparing with our predictions for the strange-baryon couplings given in Table 1, we ob- 
serve an agreement for vector-meson couplings within uncertainties. Also the convention- 
independent relative signs of T and V couplings agree. Meanwhile, the LCSR predictions 
for g\NK and gT,NK are systematically lower than the intervals for these couplings ob- 
tained in the potential model. 

3 The QGS model for meson pair production 

In the QGS model, the amplitudes of binary reactions, such as pp — )• AcAc or pp — )• 
DD, are described by planar diagrams depicted in Fig. 1. These diagrams have a dual 
interpretation. From the s-channel point of view, annihilation of the slow uu or dd 
pair from the initial proton and antiproton is followed by a creation of the cc-pair. 
The spectator quarks and antiquarks from the initial proton and antiproton coalesce 
with the created quark and antiquark to form the final state charmed hadrons. The 
intermediate state in s-channel represents a sort of a diquark-antidiquark (Fig.l a) or 
quark-antiquark (Fig.l b) string. On the other hand, in the 4-channel a virtual hadronic 
state with the quantum numbers of a charmed meson or baryon is exchanged. In the 
s ^ |t| limit, this exchange is described by the dominant Regge pole. For instance, the 
amplitude of pp — AcAc is approximated by the (degenerate) D* , D** Regge-trajectory 
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aD*{t) = a£)«(0) + a'^*i (we use the hnear approximation). The QGS-model parameters 
are obtained [21 [121 1131 [Sj using the quark-parton description of the s-channel planar 
diagram. Replacing the c-quark by s-quark in the planar diagrams of Fig. 1 we reproduce 
the QGS model for the production of strange baryons and mesons. The strange-hadron 
pair production cross section in pp collisions calculated in this model [2] agrees well with 
the experimental data. Importantly, there is a strong flavour dependence of the binary 
reactions in QGS model, encoded in the slopes and intercepts of the Regge trajectories as 
well as in the scale factors sq entering the Regge amplitudes. The relative suppression 
of the charmed hadron production corresponds, in terms of the s-channel picture, to 
a comparatively smaller probability to create a heavy quark-antiquark pair within the 
intermediate string. 

To discuss the QGS model in more detail, we first consider a relatively simple binary 
reaction involving no spins or helicities: vr^vr" — t- MM, with pseudoscalar mesons (M = 
7r^,K~^, D^) of various flavours in the final state and with isospin and/or flavour exchange 
in t-channel. The planar diagram of this process is shown in Fig. [2l At large s and 
small \t\ <C s, the scattering amplitude is written ^ in the following Regge-pole form: 

' (4) 

where s = 1 GeV^ is a universal dimensional factor and the energy dependence is 
determined by the Regge trajectory a_R(t) with the corresponding quantum numbers 
(i? = p{a2) , K*^**\ D*'-**^) . In the above, g'^^^ [t) is the residue function of the momen- 
tum transfer squared. In the QGS model [2] the F-function dependence inspired by 
Veneziano duality is adopted: 

g(-*^)(t) = C('^A^)^2r(l-aR(t)). (5) 

The coefficient C^'^*^^ is equal to the number of planar diagrams. 

The amplitude for tt~^tt~ — )• vr^vr^ is determined by the p(a2)-trajectory: 

aji{t) = ap{t) = ap(0) + a'pt . (6) 

In this case C^'^'^^ = 2, and the numerical values of the intercept ap(0) and slope a'p 
taken from [2] are presented in Table [2 

The universal parameter in QGS model can be related to the /ovrvr strong coupling 
defined as: 

(7r-(pi)7r0(p2)|p-(pi +P2)) = 9p.. e^^Hp2-piY , (7) 

where e^^^ is the polarization 4-vector of p-meson. The numerical value ^ptttt — 6.0 with 
a negligible error is then obtained from the measured [21j width: 

r(p^vrvr) = |||(p;^j3^ (8) 



where Pp^r^ = {mp/2)^\ — 4m'^/m'^ is the 3-momentum of the pions in the rest frame 
of the p. Combining the product of couplings defined in ([7|) with the p propagator 
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Figure 2: The planar diagram of tt'^-k — )• MM (q = u,s,c and M = t:^,K'^,D^). An 
additional diagram with u ^ d contributes to the 1:^11^ -production. 
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Table 2: Parameters of Regge trajectories |^ involved in the vr+vr — t- MM. 



(neglecting the width) we calculate the vr+vr — )■ vr^'vr^ scattering amplitude in a form of 
a Feynman diagram with an "elementary" p-meson exchange in t-channel. The result is 

, 2 2s + t- Ami 
"V 

At s » \t\,m1^p the above amplitude correctly reproduces the expected s"^ asymptotics, 
where J = 1 is the spin of the vector meson exchanged in i-channel. The Regge amplitude 
([1]), being analytically continued in the Mandelstam {s,t} plane to i ~ has to 
reproduce the Feynman diagram expression ([9]) at s » m^, Substituting (j5]) in ([3]) 
and expanding the F-function near t = where ap{rn?p) = 1 we obtain a pole in the 
variable t: 

which corresponds to the p-propagator pole in ([9|). Comparing the residues of the am- 
plitudes ([1]) and Q at large s, we obtain 

—^ = 2gl^^. 11 

UpS ^ 

Numerically, the above equation at the values a'p = 0.9 and (jrg/(4vr) = 2.7 adopted in 
[2] correctly reproduces the experimental value of (7p,r7r- 
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Repeating the same comparison for the binary reaction tt+tt" — >• K~^K~ with the 
strangeness-exchange in the t-channel, we replace the p trajectory by the K* trajectory 
in the Regge-pole amplitude (jl]). The corresponding parameters of GGS model are 
presented in Table [2j Note that the SU (3) violation in this model (i.e. the effect 
of a heavier s-quark) is reflected in the parameters of Regge trajectory, and also in 
the flavour-dependent normalization scale s^^, introduced in QGS approach. We then 
compare the residue function near the pole at t = rn\t with the diagram containing 
the K* propagator and the K*Ktt strong couplings. This diagram yields an expression 
similar to ([9]): 

T^l;~^''"''~\s,t) = 4^ i^s + t- 2{ml + ml) + ^"^^ " , (12) 

f mj^t — t \ ^K* J 

with the same large s asymptotic behavior. The relation analogous to (jlip yields 
gK*K-K = 4.5 for the K*^K^'k~ strong coupling. This is very close to the value ex- 
tracted from the K* — )- K-k width [21] . 

Turning to the charmed meson production in the two-pion collisions, we consider the 
amplitude (jH) with the D* Regge-trajectory: 

In the QGS approach, the flavour-dependence of the amplitude is reflected by the sub- 
stantial differences between the slope parameters of D* and p{K*) trajectories on one 
hand, and between the scale factors Sq^ and g'^'^^'^^^ on the other hand, as can be seen 
from Table El Hence as we already mentioned, there is no SU{A)fi symmetry in this 
model. Still there remains an important question if the universal value of Qq can be used 
also in the charm production amplitude. The D*Dit strong coupling is defined as0 

{TT~{pi)D'^{p2)\D*-{p^+P2))=gD'D. ef'\p2-PiY, (14) 

and the "elementary" D*-exchange diagram yields the same expression as in (jl2p . where 
K ^ D and K* — )• D* have to be replaced. Continuing the Regge amplitude to t = m^, 
and comparing with the large s limit of the L'*-exchange (jl2p . we obtain: 

-^ = 25|).D.- (15) 



Substituting the slope of the D* trajectory and using the universal value of the QGS 
model [2] we obtaiqj: gwD-w = 5.8. Interestingly, this value is close to the interval 
estimated from QCD LCSR in [22j taking into account the gluon radiative corrections 
[23j: bD*D7r]LCSR = 5.0 lb 1.75. The only existing measurement of the total D* width 
combined with the branching fraction yields a larger result: [ff_D*D7r]exp. = 8.95 it 0.15 it 
0.95 [21]. We conclude that the D* trajectory slope adopted in QGS model is consistent 
with the LCSR estimates of the D*DTr strong coupling. 



^ Note that 2gD*DTT is equal to the D*Dn couphng defined in |22) . 

^Our estimate difTers from the smaller value quoted in [IS] and based on the same model. Note that in 
this earlier paper a larger value of the slope a'^, — 0.64 was used. 
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Figure 3: Dependence of the cross sections of tt+tt^ — ?■ K^K^ (solid) and tt+tt^ — )• 
D^D^ (dashed) on piab in QGS model. 



Note that one of the inputs used to determine the D* trajectory of in QGS model is 
the Regge trajectory of J/V', taken in [2j as: 

a^(t) = -2.18 + 0.33t, (16) 

and supported by the QGS model analysis of inclusive charm production. Here one 
can mention the estimate of the intercept a^(0) obtained in [25J, where a four-point 
correlation function of heavy-quark-currents was first calculated using OPE in terms 
of loop diagram and vacuum condensates and then, via optical theorem, related to the 
photon-heavy meson scattering cross section taken in the Regge form. The comparison 
of two representations of the correlation function yields a^(0) = —(2 ~ 3), consistent 
with (fT6]l . Note that the perturbative loop approximation in [25j yields q;^(0) = 0, hence 
the estimated value of this parameter is entirely determined by nonperturbative (gluon- 
and quark-condensate) effects. 

To illustrate the QGS model for meson pair production, in Fig. [3] we present the cross 
sections of the processes discussed in this section. According to our definition of the 
scattering amplitude, the differential cross section is: 

^(tt+tt- ^ MM) = ^- -AM^ , (17) 

dt^ ' 167rs(s-4m2) ' ^ ' 

where we substitute the Regge-pole amplitude ^ for M = K,D and integrate over \t\ 
from the kinematically allowed minimal |to| to |to| + A. Here we choose for definiteness 
A = 0.6 GeV^, so that \t\ remains much smaller than s and hence the Regge-pole descrip- 
tion is applicable. Hereafter we refrain from predicting total cross sections, because at 
large |t| ~ s the behavior of the scattering amplitudes is governed by mechanisms other 
than a simple Regge-pole model. Still, differential cross sections rapidly decrease with 
\t\, hence, our results for the integrated cross sections provide an order-of-magnitude 
estimate also for the total cross sections. As we can see from Fig. [31 the charmed meson 
production cross section is strongly suppressed with respect to the strange-meson one. 
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4 pp- production of charmed and strange hadrons 

The QGS model of the charmed baryon-pair production in pp colhsion is described by 
the planar diagram in Fig. la. This amplitude, similar to vrvr — >• DD, is approximated 
by the D* Regge-trajectory. The amplitude pp — )• Ac Ac presented in [2] has a helicity- 
averaged form 

' ^^^^ 

where the residue function g^^^^^t) = CP^'=5gr(l — aD*it)) contains the same universal 
coupling (7q . Importantly, the scale factor Sq^" obtained following [2] (see Table [3]) is 
not equal to Sq^, reflecting the difference between baryon and meson production in this 
model. 

Here we consider a more elaborated version of the QGS model for this process with 
the helicity amplitudes (see App. A). The differential cross section has the following 
expression 



^{pp^AcAc ^ 



dt 327rs(s — 4m|) 

+2\H{++, + \H{—, + |F(-+, + |F(+-, 



//(++, + 2|F(+-, 

(19) 



where in the helicity amplitudes i/(AiA2; A3A4) the notation Ai^2; (^^3,4) denotes the he- 
licities of the proton and antiproton (Ac and Ac), respectively. The s,t dependence of 
the amplitudes is not shown for brevity. At fixed s, the region of the momentum transfer 
squared t is given by ti < t < to; where: 



*o(i) = + - 2 + (-) 2 V " '^'^P^^^ ~ ^""'c) ■ (20) 

We assume that each helicity amplitude has the Regge form (jlSp . The residue functions 
are fixed by continuing these Regge amplitudes to the point t = m^, and matching them 
to the helicity amplitudes of pp — )• AcAc with an "elementary" D* exchange. The latter 
are given in App. A. (see eq.(|3ip) and contain two independent strong AcD*N couplings 
defined in ([T|), so that the coupling g\^^D* {oji^ND*) enters the helicity-nonflip (-flip) 
amplitudes. We arrive at the following expression for the cross section: 



(211 



where the function 



Rd' (s, t) = a'B.r(l - an* (t)) s i—j (22) 
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is determined by the Regge-pole parameters. As opposed to vrvr — )• MM, where we 
used the original QGS model with the universal normalization parameter Qq, there is 
now a more subtle substructure of the Regge amplitudes with the strong couplings 
determining the helicity-flip and helicity-nonflip amplitudes. Furthermore, we modified 
the above cross section with respect to ()19p by multiplying it with the so called absorption 
factor (^fr'P~^^cAc) ^j^jg]^ jg included following [2j. This factor derived in App. B takes 
into account the initial- and final-state rescattering of the baryons and antibaryons and 
suppresses the cross sections. 

The related processes of charmed baryon production: pp — ScAc and pp Sc^c, 
have a similar description in QGS model, in particular, they are also dominated by 
the same D* Regge-pole exchange. Their cross sections depend on the combinations 
of couplings Qj^ j^jjt: and 5^' . The corresponding expressions in terms of helicity 
amplitudes have minor difi^erences with respect to ()19p which we will not discuss here 
for brevity. The numerical analysis yields a substantial suppression of the Sc production 
versus Ac production, due to the difference in the strong couplings inferred from LCSR. 
This suppression will be discussed below in more detail. 

The charmed- meson production, pp — )• DD, is described in QGS model by the planar 
diagram depicted in Fig. lb. In this case the t-channel exchange involves Ac and Sc 
Regge-trajectories. Their parameters presented in Table [3] are assumed equal. However 
according to our predictions, the strong couplings of these baryons to mesons and nucle- 
ons quite differ from each other, hence there is a significant difference between the cross 
sections of charged and neutral charmed meson-pair production. Indeed, in the planar 
diagram model, the process pp — )• D^D^ can only be mediated by the S^"*" exchange in 
t-channel, whereas in pp — J- D^D^ both trajectories Ac and enter the amplitude. Note 
that this is a characteristic feature of the planar diagram mechanism. Inelastic scatter- 
ing in the final state (D^D^ — )■ D~^D~) due to nonplaner diagrams can enhance D~^D~ 
production cross section. Moreover, in a model where these processes are mediated by 
intermediate charmonium states in s-channel, pp — )• {cc} — )• DD, there is no correlation 
between the flavours of initial and final hadrons, so that both charged and neutral D 
mesons are produced with equal rates. Such a model may indeed work as an additional 
mechanism for charmed meson-pair production slightly above the threshold, (see e.g., 
|10j ) but the resulting cross section is much smaller than the one generated by i-channel 
exchanges. Let as also mention that, according to the model [9] based on the baryon- 
antibaryon potential, the initial state inelastic interaction could significantly enhance 
the D^D~ production cross section in the near-threshold region. Therefore, the charged 
charmed meson cross section can serve as a useful check of different charm-production 
models. 

The decomposition in the helicity amplitudes in pp — )• DD is simpler than for the 
baryon-pair production because only the helicities of the initial proton and antiproton 
are involved. For the D^D^ production we follow the same method of matching the 
Regge-pole amplitude to the "elementary" Ac-exchange at t = m^^ and obtain the cross 
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Table 3: Parameters of the Regge trajectories determining the pp amplitudes of charmed 
and strange hadron-pair production in QGS model J^. 



section: 
where 

i?A,(s,t) = a:,^r(--aA.(t)) f-^j , (24) 

and the Sc exchange contribution is neglected due to much smaller couplings. The 
absorption factor C^^^^ ^ ^ in the above cross section has a form similar to the one in 

m- 

Turning to the numerical analysis of the cross sections we notice that LCSR predic- 
tions for strong couplings have a typical error of ~ 50%, hence their fourth powers in 
the cross sections introduce large uncertainties. This mainly concerns the ^(^-couplings. 
The ratios j are predicted from LCSR with much smaller uncertainties and more- 
over, the helicity-flip contributions to the cross sections proportional to 5^-couplings are 
kinematically suppressed at small t. 

In order to decrease the uncertainty of the predicted cross sections for charmed hadrons 
we consider also the strange hadron pair-production in pp collisions. Extending the 
(modified) QGS model to these processes, allows us to test it, because in this case some 
(albeit, quite old) experimental data are available. Moreover, we use the fact that the 
ratios of the strong couplings of charmed and strange hadrons given in Table [T] have 
comparatively smaller uncertainties, than the individual couplings, because the same 
nucleon DA's are used in the sum rules in both cases of charmed and strange hadrons. 
Hence, we can constrain the couplings of strange hadrons by fitting the model cross 
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Figure 4: Differential cross sections of pp — t- AA and pp — )• SA at piab = 6 GeV, and 
pp — )■ K~^K~ at piab = 4 GeV. The data points are from The solid 

curves are given by QGS model, with the ratio of tensor to vector strong cou- 
plings from LCSR (dashed curves indicate the uncertainties) and the vector 
strong couplings fitted to the measured cross-section normalization. 



sections to experimental data and then use the calculated ratios of the couplings to 
reproduce the charm cross sections with smaller uncertainties. 

Let us start with the process pp — )• AA. Its cross section in the QGS model has the 
same form as ([2T]) . with the strong couplings ff^JViC*' ^* Regge-trajectory and the 
absorption factor c^^^^^^ . We first calculate the differential cross section ^ {pp — )• AA) , 
without taking into account the absorption factor, and fit the slope of the t-dependence to 
the exponential form exp{—Lji(s)\t\). In this cross section we use the ratio of tensor and 
vector strong couplings, g'j^jqx* / d^NK* ' obtained from LCSR (see Tabled]) . The slope Lr 

is then used to calculate the absorption factor c^^^^^^ as explained in App. B. Note that 
the overall normalization of the cross section depending on the vector strong coupling 
g^NK* does not play role in this determination. On the other hand, due to difference 
of the slopes for the Regge amplitudes of strange and charmed hadron production 
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the resulting absorption factor turns out to be almost twice larger for pp — t- AA than for 
pp — )• Ac Ac at small t (see App. B), in accordance with the estimates in 

After inserting the calculated absorption factor in the differential cross section of 
pp — )■ AA, in FigU] we compare the latter with the data points [26] at piab = 6 GeV 
and at small t where we expect the QGS model to work. Note that not only the Regge 
amplitude itself but also the absorption factor contribute to t-dependence, making it 
steeper. As can be seen from this figure, the agreement of the shape of the differential 
cross sections with the data is not very good, which can possibly be traced back to slightly 
oversimplified model of t dependence for this particular (not yet sufficiently large) energy 
in our model. Still in the integrated cross section which is our main interest here, we 
expect that the imperfection of the shape does not play an important role. 

As a next step, we fit the overall normalization of the cross section to the data and 
obtain the interval of the vector strong coupling 

I^AATif*! = S.StoJ, (25) 

which is within the broader interval of the LCSR prediction given in Table [H After 
that, we combine the above estimate with the calculated ratio gX^j^D* /9\nk* (^^^ 
bled]) and estimate the vector coupling for the charm case |(7a^7vd*I ~ ^-^^le) again in 
agreement with the interval of LCSR prediction. We use the above "rescaled" interval 
for g\^^£i* in obtaining the charmed baryon cross section ()2ip . thereby decreasing the 
resulting uncertainty. Note that the ratio Oa^nd* /dXcND* again taken from the LCSR 
prediction. The cross section of Ac pair production we are interested in is then calcu- 
lated in two steps: first we fix the exponential slope Lji in order to obtain the absorption 
factor and second, include this factor in the cross section. To estimate the cross sections 
of pp — >• Tic Ac, ScSc and pp — >• we use a similar procedure employing the available 

data on strange hadron pair production (see Fig. In particular, fitting of the corre- 
sponding strong couplings yields l^sAfit*! ~ ^-^toi and \gANK\ = 13.9^07- The LCSR 
predictions for these couplings given in Table [T] are only marginally consistent with the 
above intervals. Note that the predictions of the potential model [20j with the scatteing 
potentials obeying a (slightly broken) SU{3) fi symmetry, for the same couplings are in 
a better agreement with the fitted values. 

Differential cross sections of pp — t- AcAc, Sc-Ac, ScSc and pp — t- DD are displayed in 
Fig. [5] as a function of t at piab = 15 GeV. As expected, their slope is much smaller 
than in the case of strange hadron production. The integrated cross section a{tQ,A) of 
charmed baryon or meson pair-production is defined as the integral of the differential 
cross section over the region of small momentum transfers: max{ti,to — ^} < t < to, 
where we adopt A = 0.6 GeV^. These cross sections plotted as a function of piab in the 
region accessible to PANDA are presented in Fig. [6l The summary of our results for 
the cross sections is also displayed in Table HI 

Let us emphasize that the uncertainties of the predicted cross sections are still quite 
large, even after we narrowed them using the strange hadron production data. Note 
that we only quote the uncertainties stemming from the LCSR estimates of the strong 
couplings. The QGS model itself has "systematical" uncertainties, which is difficult 
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Figure 5: Differential cross sections of pp — )• AcAc, and pp — )• DD at piab = 15 GeV 
calculated in QGS model. The dashed lines indicate the uncertainties caused 
by LCSR estimates of strong couplings. 



to assess quantitatively, as it is the case for any phenomenofogical hadronic model not 
directly related to QCD. The predictive power of the model concerns mostly the ratios of 
cross sections, where the "intrinsic" uncertainties of the method to a large extent cancel. 
One important prediction concerns the suppression of Sc- with respect to Ac-production 
cross section. 

This suppression is more significant than predicted in [2] where simple relations based 
on the nonrelativistic quark-diquark model for these reactions are used 



a{pp Ac Ac) _ ajpp ScAc) 
a{pp ScAc) a{pp Sc^c) 

Our predictions for these ratios at piab = 15 GeV are: 



3 . (26) 



ajpp Ac Ac) _ +1.0 (t(pp ^ ScAc) _ +o.9 

a{pp ^ ScAc) - a{pp ^ ScSc) - ■ ^'^^ 

Due to the suppression of Sc couplings versus Ac couplings, also the D^D^ production 
cross section is expected to be significantly larger than the D~^D~ one. It will be very 
interesting to test all these predictions experimentally. 



5 Conclusion 

In this paper we bring together the strong couplings of charmed and strange baryons, 
both predicted within one and the same QCD based method of LCSR. We have demon- 
strated that it is possible to avoid SU (4) fi approximation. The relations between cou- 
plings are nontrivial because they stem from the nonperturbative dynamics which is 
quite different for heavy and light (also strange) quarks. The LCSR predictions for 
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channel 






pp AcAc 


130 (30 ^ 470) 


60(15^210) 


pp ScAc 


24(5.0^ 140) 


12(2.0^70) 


pp ScSc 


5.0(1.0 ^45) 


3.0(0.4^24) 


pp D^D^ 


52 (13 ^ 200) 


20(5.0^75) 


pp — ^ D^D^ 


< 0.01 


< 0.01 



Table 4: Differential and integrated cross sections with A = 0.6 GeV^ for charmed 
hadron production at piat = 15GeV. 

strong couplings can be significantly improved in future by calculating radiative gluon 
corrections and taking into account soft gluon components of the nucleon DA's. 

The main task of this paper is to estimate the charm production cross sections in pp 
collisions. For that purpose we have selected the most (in our opinion) "QCD-friendly" 
model of hadronic reactions, that is, the Kaidalov's QGS model. This approach has 
revealed itself as a very useful tool for hadronic reactions with different flavours, also for 
inclusive production of hadrons. In this paper we used a more detailed description of 
binary processes with baryons in terms of helicity amplitudes and employed the strong 
couplings of initial protons and final charmed baryons (mesons) with the intermediate 
charmed mesons (baryons) calculated from LCSR. 

Strictly speaking, the QGS model is applicable only at sufficiently large energies, 
beyond the upper limit of the PANDA energy region. Hence the cross sections calculated 
here can only be considered as an order of magnitude estimates, also because the model 
is only valid at small momentum transfers and the absorption factor is only taken in the 
first approximation. Still the relations between cross sections are less influenced by the 
uncertainties and are almost independent of the absorption factors. In future, the model 
adopted in this paper can be developed further, taking into account of the subleading 
Regge trajectories and a more elaborated absorption ansatz. 

Finally, turning to the comparison of our results with the charm-production estimates 
in the literature, we observe that our prediction for the dominant AcAc production is 
(within estimated uncertainties) consistent with the one obtained in the original QGS 
model [2]: a{pp — )• AcAc) ~ 100 nb, at piab = 15 GeV, whereas the predictions for the 
ratios of cross sections obtained here and in [2] differ. For example, we do not exclude 
a larger charmed meson cross section than a{pp — t- D^D^) ~ 5 nb predicted in [2]. 
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Figure 6: The integrated cross sections cr{to, A) of charmed baryon and meson pair pro- 
duction in pp collisions in QGS model. The dashed lines indicate the uncer- 
tainties introduced by the strong couplings obtained from LCSR. 



A model of exclusive charm production cross sections based on the QGS and Regge- 
poles can be found in [3], where the SU{^)fi symmetry is used and a different form of 
the cross section is adopted, adding a t-dependent dipole "residual factor". Numerically, 
our predicted intervals for the differential cross sections at to turn out to be larger than 
the ones in [3J. 

The other models in the literature are based on radically different approaches. E.g., in 
[8^ a hadronic baryon-antibaryon potential derived from the coupled channel approach is 
used, predicting the cross section of AcAc production up to a few //6 near the threshold, 
i.e., much larger than obtained here. On the opposite side are the typically smaller 
cross sections obtained from perturbative approaches, such as the inclusive charm pro- 
duction estimate in the parton model [7] and the approach [6j to pp — t- AcAc employing 
distribution amplitudes of initial and final baryons. 

Concluding, this paper contains an attempt to apply QCD predictions for hadronic 
strong couplings to the models of exclusive hadronic reactions. Our estimates for charm 
production cross sections contain rather large uncertainties. Still even the lower limit of 
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the cross sections predicted here allows one to expect an appreciable number of charmed 
baryons and mesons produced at PANDA. We look forward to other applications of 
the strong couplings presented in this paper and their future improvements. 
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Appendix A: Helicity amplitudes 

The helicity amplitudes of pp — )• AcAc scattering via D* meson exchange are obtained 
from the initial invariant scattering amplitude 



H{XiX2; A3A4) 



1 



m 



D* 

%(P2, A2 



T 



rriA^ + TUN 



un{pi, Ai 



T 

V J* , ■ SAcND* 



*P T 



vaSpa,X4), (28) 



where the baryon bisponors are distinguished by their indices so that (pi,Ai), (p2,A2), 
(PSj-^s) and (p4,A4) are the four-momenta and helicities of the proton, antiproton, Ac 
and Ac respectively; is the polarization vector of the virtual D* meson. Generally, 
there are 16 different helicity amplitudes for pp — )• AcAc process, only six of them are 
independent due to symmetries PH] . 

In the cm. frame of proton-antiproton pair we choose the x, z plane for the process, 
and the 3-momentum of proton in the z direction, so that the 3-momentum of Ac has 
the angular coordinates {6, 99 = 0). Then the kinematics is as follows: 



1 



1 



Pi = 2 0, 0, - 4m2,) , p2 = -{Vs, 0, 0,-^Js- Am%) , 



P3 



{\/s , \/ s — 4m\ sm6,0, a/ s — 4m\ cos 9), 



P4 = 2 (v^' ~ V * ~ 4"^Ac sin 6*, 0, — ys — 4m^^ cos ( 



Summing over the polarization of D* meson: 

^ eM(g^A)e-((/,A) = -5^'' + 



A=l, 2,3,4 



m 



2 ' 

D' 



(29) 



(30) 



and substituting explicitly the bispinors with various helicities in the chosen frame we 
obtain rather bulky expressions of helicity amplitudes, which however greatly simplify in 
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the limit of large s where we compare them with the Regge amplitudes. The six helicity 
amplitudes are: 



H{+- 
H{++ 
H{— 



H{-+ 
H{+- 



++ 



+ 



++ 
-+ 

-+ 



t — m 



D* 



t — m 



D* 
s 



ND* 



'^\^9a_^nd*9AcND* 

'^V-tdA^ND'dA^ND* 



t — m 



D* 



t — m 



D* 



t - m|,. (ttt-a, + mAr)2 ' 



(31) 



It is clear that in the large s limit only three helicity amplitudes are independent. The 
amplitudes can be related to each other through the following relations: 



H{—Xi — A2; —A3 — A4) 
ff(A2Ai; A4A3) 



-l)^i-^2-^^^+^*F(AiA2;A3A4), 
-l)^i-^^-^^+^^//(AiA2;A3A4), 



(32) 



following from the parity and charge-conjugation invariance. 



Appendix B: Absorption factor 

Here we derive the absorption factor C^(s, t), multiplying the cross section. The absorp- 
tion is generated by the (quasi) elastic rescattering of the initial proton and antiproton 
as well as of the final hadron pair, both are approximated by the pomeron exchange [2]. 
Here we make a simplifying assumption that the elastic rescattering amplitudes domi- 
nate, they do not change the helicities and are the same in the initial and final states, 
independent of the flavour content of the latter. 

Consider a process pp — t- BB with generic B hadrons in the final-state. The amplitude 
in QGS model, having the form (jlSp has predominantly exponential behavior at small 
t, the main source of it is the Regge-pole factor (s/sq)"^*^- Therefore, the pp — )• BB 
amplitude can be cast in the exponential form 

TR{s,t) = /i?,(s)exp • (33) 

We then switch to the impact parameter representation, where the 2-dimensional vector 
b is conjugate to the transverse momentum transfer q±: 

Tnis, h) = j Tr{s, t) exp(ig-l • 5)^ , (34) 
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and at high energies t = ~ "l^l^- The angular integration yields: 

oo 

Tr{s, b) = ^j d\t\Jo (^^At\b) Tr{s, -\t\) , 



(35) 



where Jq is the Bessel function and b = \b\. Substituting the exponential representation 
()33p in the above integral, and integrating over t we obtain 



TR{s,b) 



Lr 



■ exp 



(36) 



The rescattering in the initial and final state is dominated by the pomeron amplitude 
Tp(s,t) which is predominantly imaginary and has an exponential form in the momen- 
tum transfer Tp{s, t) = Tp{s, 0) exp(— Lp|t|/2). The forward-scattering amplitude is ex- 
pressed via total pp cross section using the optical theorem: ImTp(s, 0) = 2p*^a^°p{s), 
where p* is the 3-momentum in the cm. system of the pp collision. Hence one obtains 
for the pomeron-mediated elastic rescattering amplitude: 



Tp{s,t) = 2ip*^sa';^{s)eM-Lp\t\/2). 

Employing the same Fourier-transformation to the impact parameter space as in 
it is easy to get the impact parameter representation for this amplitude: 



(37) 



Tp{s,b) 



2ip*^ 



Lp 



exp 



2Lp 



(38) 



The absorption contribution added to the initial Reggeon amplitude in the b space 
yields: 

Tpis,by 



Tis,b) = TR{s,b) 



1 + i 



S-Kp*y/S 



TR{s,b) 



1-X{s,b) 



where 



Xis,b) 



a. 



is) 



■ exp 



(39) 



(40) 



and the normalization factor multiplying Tp corresponds to the convention of impact 
parameter representation adopted in p]. Substituting ()36p in ()39p and performing the 
inverse Fourier transformation to the t-dependent amplitude we finally obtain: 



T{s,t) = TR{s,t) 



1 



tot 

pp 



is) 



Ll\t\ 



(41) 



This expression takes into account absorption in the amplitude in the first approxima- 
tion. To obtain C^(s, t) one simply has to square the expression in brackets multiplying 
Tp{s,t) in the above. At high energies this effect should be resummed (exponentiated), 
however at small i's we are considering here the resummation effects, as we checked 
numerically, are small, so that the first approximation for Ca is sufficient. 



20 



For the numerical evaluation of the absorption factor in (j4ip . the data on cr^{s) in a 
parameterized form are taken from [21j, so that crj^p{s) changes from 52.5 mb to 47.9 mb 
in the interval piab = 10 GeV to 20 GeV. The slope of the pomeron mediated elastic 
pp scattering is taken from [29j, e.g., Lp = 12.1 GeV~^ at piab = 15 GeV. Finally, 
the slopes of Regge-pole amplitudes fitted to the exponential form (|33p are (in units 
GeV-2): Lr = 2.6 {piat = 6 GeV) and _2.5 {piat = 4 GeV) for pp AA, SA and 
K^K~, respectively. For pp — )• Ac Ac, ScAc, ^c^c, DD, the corresponding slopes are 
Lji = 0.6, 0.4, 0.2, 1.2 {piab = 15 GeV), respectively. For numerical illustration, we 
quote the absorption factors calculated at the same eiiergy piab = 15 GeV and at t = to 
for strange and charmed baryon production: c^p^^^^ = 0.09 and = 0.04. 

These values are in agreement with (t-averaged) values presented in [2j and indicate a 
strong absorption effect on one hand and a large difference between this effect for strange 
and charmed baryons. 
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